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Abstract Charged systems with partially annealed charge disorder are investigated using 
field-theoretic and replica methods. Charge disorder is assumed to be confined to macroion 
surfaces surrounded by a cloud of mobile neutralizing counterions in an aqueous solvent. 
A general formalism is developed by assuming that the disorder is partially annealed (with 
purely annealed and purely quenched disorder included as special cases), i.e., we assume 
in general that the disorder undergoes a slow dynamics relative to fast-relaxing counterions 
making it possible thus to study the stationary-state properties of the system using methods 
similar to those available in equilibrium statistical mechanics. By focusing on the specific 
case of two planar surfaces of equal mean surface charge and disorder variance, it is shown 
that partial annealing of the quenched disorder leads to renormalization of the mean surface 
charge density and thus a reduction of the inter-plate repulsion on the mean-field or weak- 
coupling level. In the strong-coupling limit, charge disorder induces a long-range attraction 
resulting in a continuous disorder-driven collapse transition for the two surfaces as the dis- 
order variance exceeds a threshold value. Disorder annealing further enhances the attraction 
and, in the limit of low screening, leads to a global attractive instability in the system. 
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1 Introduction 

Interaction of charged macromolecules (macroions) is essential for soft and biological ma- 
terials in order to maintain their complex structure and distinct functioning. In many cases, 
charge patterns along macromolecular surfaces are inhomogeneous and exhibit a highly dis- 
ordered spatial distribution. DNA microarrays [1], surfactant-coated surfaces (2), random 
polyelectrolytes and polyampholytes |3 1 present examples of such disordered charge dis- 
tributions. The charge pattern can be either set and quenched in the process of assembly 
of these surfaces, or can exhibit various degrees of annealing when interacting with other 
macromolecules in aqueous solutions. Disorder annealing in charged systems may result 
from different sources; e.g., finite mobility and mixing of charged units (lipids and proteins) 
in lipid membranes [4], conformational rearrangement of DNA chains in DNA microar- 
rays 1 1] and charge regulation of contact surfaces bearing weak acidic groups in aqueous 
solutions 1 5 1, to name a few, all lead to annealing effects. In reality, one may deal with a 
more complex situation where the surface charge pattern displays an intermediate character 
[2] and thus may neither be considered as purely quenched (i.e., with fixed random spatial 
distribution) nor as purely annealed (i.e., thermally equilibrated with the bulk solution). 

Charge disorder appears to produce electrostatic features that are remarkably different 
from those found in non-disordered systems. Mounting experimental evidence shows that 
like-charged phospholipid membranes and fluid vesicles, which primarily contain mobile 
surface charges, may undergo aggregation and fusion in the presence of multivalent cations 
131. A similar behavior is observed with negatively charged mica surfaces exposed to a so- 
lution of positively charged surfactants (2J; here formation of a random mosaic of surfactant 
patches on apposing surfaces (after the surfactant is adsorbed from the bathing solution onto 
the surfaces) leads to a long-range attraction and thus a spontaneous jump to a collapsed 
state. Although this behavior is akin to the transition to the primary minimum within the 
standard DLVO theory of weakly charged systems [6 7], the attractive forces at work here 
exceed the universal van-der-Waals forces (HJ incorporated in the DLVO theory by a few 
orders of magnitude |2|. 

In fact, the emergence of an instability is not captured by the standard theories of charged 
systems that incorporate static, non-disordered charge distributions for macromolecular sur- 
faces. These theories cover both mean-field [6.7, 9] or weak-coupling limit (including the 
Gaussian-fluctuations correction around the mean-field solution flO|) as well as the strong- 
coupling (SC) limit [11 12, 13, 14], where the central theme is the absence or emergence 
of electrostatic correlations induced by neutralizing counterions in the system that give rise 
to attractive interactions between like-charged objects. Both uniform |6 7,9 12] as well as 
modulated [15] charge distributions have been considered in this context. In the mean-field 
regime, like-charged objects always repel. While the opposite limit of strong coupling (real- 
ized, e.g., with high valency counterions, highly charged macroions, low medium dielectric 
constant or low temperature I11II12I ) is dominated by correlation-induced attractive forces 
that can bring two apposing like-charged surfaces to very small separation distances. At 
small separations, however, a universal repulsion due to the confinement entropy of inter- 
vening counterions sets in and stabilizes the system in a closely packed bound state with 
a finite surface-surface separation 11111121 . This is true even for surfaces of opposite (un- 
equal) uniform charge distribution [16]. Therefore, other mechanisms have to be at work 
that would lead to attractions strong enough to counteract such repulsive forces and lead to 
collapse or instability in a system of charged macroions. One such mechanism we propose is 
the disorder of the charge distribution along macromolecular surfaces that turns out to be as 
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significant as the counterionic correlations and could provide a new paradigm in the theory 
of charged soft matter. 

Previous studies of charge disorder on macromolecular surfaces have investigated both 
types of quenched 1 3 ,17, 18 19] and annealed [3, 19,20,21 ,22,23] disorder (including, specif- 
ically, the classical work on charge-regulating surfaces [5]). They mainly deal with situa- 
tions where the system is in equilibrium and, on the question of electrostatic interaction 1211 
22.23 1, focus primarily on the weak-coupling regime, where disordered surfaces of equal 
mean charge always repel and no collapse or instability arises [24]. A systematic analysis of 
quenched disorder effects is presented in the previous works of two of the present authors 
11711 181 . where it was shown that not only can electrostatic interactions between like-charged 
objects turn from repulsive to attractive due to counterionic correlations, but also the disor- 
der of the surface charge itself can give rise to an additive long-range attraction. This is 
most clearly demonstrated by attraction induced between disordered surfaces of zero mean 
charge but with a finite variance of the disordered charge distribution [18]. In the SC limit, 
the quenched disorder-induced attraction may be so strong that it can dominate the entropic 
repulsion at small separations and continuously shrink the SC surface-surface bound state 
1 1 1 II 121 upon increasing the quenched disorder variance, predicting thus a continuous col- 
lapse transition between a stable and a collapsed phase beyond a threshold disorder variance 
1171 . Since the experimental situation may be more complex [2], and it might not allow for 
straightforward differentiation of the charge pattern into purely quenched and purely an- 
nealed cases, we next set ourselves to explore possible effects from partial annealing of the 
surface charge. If there is a fingerprint of the partially annealed surface charge disorder on 
the nature and magnitude of surface interactions, this would help in assessing whether the 
experimentally observed interactions can be interpreted in terms of disorder-induced inter- 
actions or not. This is the motivation with which we venture on this exploration. 

In this paper, we present a general formalism for charged systems with partially an- 
nealed disorder by invoking field- theoretic and replica methods. We then focus on the case 
of two interacting planar charged surfaces as a model system and examine explicitly the 
effects of disorder on the inter-surface interaction in this system. Partially annealed disorder 
in general arises when a coupled motion of slow and fast variables (corresponding here to 
surface charges and counterions, respectively) is present in the system. It represents a non- 
equilibrium situation, whose investigation requires suitable methods. The previously stud- 
ied cases of static, non-disordered surface charge distribution I11U12I as well as quenched 
1171 and annealed surface charge distribution follow as special cases from our formalism. 
In the SC limit, we find that the system of two like-charged planar surfaces with neutral- 
izing counterions becomes globally unstable upon annealing the quenched surface charge 
and collapses into contact regardless of other system parameters due to strong attractive 
forces from the annealing effects. Hence, the quenched phase behavior is not stable against 
small annealing perturbations and is dramatically changed. However, stability may be re- 
stored in this system by adding a finite amount of added salt that screens out the long-range 
Coulomb interactions. In this case, we recover the continuous collapse transition between 
a stabilized closely packed bound state of the two surfaces and a collapsed state where the 
surfaces are in contact. This is qualitatively similar to the purely quenched case | "17"]. How- 
ever, the partially annealed bound state shows a significantly larger attraction and a smaller 
optimal surface separation as compared to the quenched case. In other words, allowing for 
rearrangements of the macroion charges leads to configurations of lower free energy. Since 
the present formalism is quite general, we shall also study the mean-field limit, where (in 
contrast to the quenched case where no disorder effects are found |17,19,24|) the disor- 
der annealing appears to suppress the mean-field repulsion significantly by renormalizing 
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Fig. 1 Schematic view of a system of macroions with partially annealed disordered charge distribution p (r) 
and ^-valency counterions at bulk temperature T. Surface charges may exhibit a different effective tempera- 
ture T' due to their disordered nature and slow dynamics relative to the fast-relaxing counterions. As a model 
system, we study two apposed planar macroion surfaces with disordered surface charge distributions (speci- 
fied in the text) located at z = ±a at the separation distance d = 2a. We neglect the dielectric discontinuity at 
the boundaries or the presence of added salt in the system 1 18 1 (see also Refs. [25, 26 1271 ). 



the surface charge to smaller values. Hence, besides the previously established mechanisms 
of counterion-induced [12 13 14 1 and quenched disorder-induced II 1 "711 181 correlations, we 
find that the annealing of macroion charges provides another mechanism enhancing the like- 
charge attraction. 

The organization of the paper is as follows: We start with the general formalism that 
allows us to define and to deal with the partially annealed disorder in terms of an "effective 
partition function" that is obtained in the form of a functional integral over a fluctuating local 
electric potential field. The structure of this field theory is too complicated to allow for a 
general solution. We thus derive asymptotic solutions in the mean-field limit (corresponding 
to the Poisson-Boltzmann theory of the classical DLVO framework) as well as the strong- 
coupling limit via an application of the replica formalism. We finally evaluate and analyze 
the inter-surface interactions for the specific case of planar charged surfaces in both limits 
and compare them. We conclude by positioning our results in the growing framework of the 
weak-strong coupling formalism for charged macromolecular interactions. 



2 General formalism 

Let us consider a system of fixed macroions with disordered charge distribution, p(r), im- 
mersed in an aqueous medium of dielectric constant, e, along with their point-like neu- 
tralizing counterions of valency q. In what follows, we shall develop our formalism for an 
arbitrary ensemble of fixed macroions but for explicit calculations, we shall delimit our- 
selves to a model system of two apposed charged planar surfaces with p (r) representing the 
charge distribution along both planar surfaces (see Fig.QJ. In the quenched limit, p(r) is as- 
sumed to be static and only counterions are subject to thermal fluctuations. In the annealed 
limit, both counterions and macroion charges are subject to fluctuations of comparable time 
scales {i.e., T c j ~ T s respectively) and thus mutually equilibrate. The intermediate situation 
of partially annealed disorder by definition occurs when there is a macroscopic separation 
of time scales between the so-called fast and slow variables as frequently observed in glassy 
systems I28II29I . 

In the present context, counterions comprise the fast variables as they are dispersed 
and freely fluctuate in the bulk. Macroion surface charges are assumed to constitute the 
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slow variables with the time scale T s s> T c j due to their disordered nature (as, e.g., they are 
confined typically within closely packed or quasi-two-dimensional disordered regions such 
as in lipid bilayers |4], surfactant-coated surfaces or surface hemimicelles 030)). 

Under these conditions, the mutual equilibration of fast and slow variables is hindered. 
Counterions rapidly attain their equilibrium at bulk temperature T and thus, because of 
the wide time-scale gap, their equilibrium free energy acts as a driving force pushing the 
slow dynamics of the surface charges to reach a non-equilibrium stationary state at long 
times. This scheme, known generally as the adiabatic elimination of fast variables 1 3 132 
[33), is investigated in a growing number of works, for instance, in the context of far-from- 
equilibrium stationary states and thermodynamics of two-temperature systems [34, 35)- It 
has been applied in particular to study spin glasses with partially annealed disorder of the 
spin-spin coupling strength 1 28 29 , 36 ] . It has been shown in general that the stationary state 
of such systems may be described by a Boltzmann-type probability distribution featuring the 
temperature of fast degrees of freedom T as well as an effective temperature T' associated 
with the disorder. This peculiar two-temperature representation clearly reflects the intrinsi- 
cally non-equilibrium nature of partial annealing. Obviously, the equilibrium free energy of 
fast variables {i.e., counterions in our case) will show up explicitly in the aforementioned 
probability distribution (here we shall not consider the relaxational dynamics of the system 
and focus only on stationary-state properties). 

In general, one may thus define an effective "partition function", 3f, in analogy with 
the equilibrium partition function that greatly facilitates the analysis of the system far from 
equilibrium [28 ,29 . 35 36 1. This procedure is discussed in Appendix[A]by adopting a simple 
dynamical model for a charged system with surface charge disorder and by identifying its 
stationary-state probability distribution. We thus find 

JT = J $>pexp(-f3'W[p]), (1) 

where /3' = j^jr and the density functional W[p] can be cast into the form 

P'np] = ~/ drg-^r) [p(r) -p (r)] 2 -n In 2f ci [p], (2) 

where 

T 

n = -. (3) 

Equation |2]l includes statistics of both counterions and the disordered charges on macroion 
surfaces. The first term is the contribution of the disorder. It can be interpreted as a general 
effective disorder potential expanded to the second order around a typical value po (Ap- 
pendix[A]i, which is always possible if one interprets g(r) as playing the role of an effective 
disorder "compressibility". This expansion leads to the standard Gaussian disorder weight 
with the mean value po(r) and variance g(r) that can be handled most conveniently by 
replica techniques 1 28.29]. Namely, 

sr = J 9p&\p]fa\p]Y = lfa\p]\ H} j, (4) 

where double-brackets denote the average {(•■•)) = / @p£?[p] (■■■) with respect to the 
Gaussian probability distribution 

9>\p\ = C exp (- l - fdrg- 1 (r) [p(r) - p (r)A (5) 
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with C being a normalization factor. 

The second term in Eq. 101 is the equilibrium free energy of a system of counterions 
at & fixed realization of disordered macroion charge, p = p(r). It follows by integrating 
over the counterionic degrees of freedom equilibrated at temperature T . In grand-canonical 
ensemble, the fixed-p partition function, 3f c {[p], can be cast into a form of a functional 
integral as 11011121 

% 

where (j) (r) is the fluctuating electrostatic potential field, jS = -Aj and 



^ci[p] = / -ar 



(6) 



dr 



ego 
2 



V0) + ip<p-lk B TQ(r)e 



-iPqe (j> 



(7) 



is the effective Hamiltonian of the system comprising Coulomb interaction v(x) = (4^eeo|x|) _1 
between all charged units (the first two terms) as well as the entropy of counterions (the last 
term). Here A is the fugacity, 3f v = \J detjSv(r,r'), and X2(r) is a geometry function that 
specifies the free volume available to counterions, i.e., the space between the two apposed 
planar surfaces in the model system in Fig.Q] 

The partition function (T} can be evaluated by using the replica trick II28U29I , i.e., by 
taking n an integral number and then standardly extending the results to real axis (for any 
real value of n = jS'/jS) by analytical continuation. Thus by using Eq. ((SJ and averaging over 
p (r) , we arrive at the disorder-averaged expression 



where the n-replica effective Hamiltonian reads 



a.b 



^[{<M] = ~E / drdr> a (r)S a6 (r,r , )fe(r') + 



(8) 



(9) 



+ £ / dr 



i J Spo(r)^(r)-AX2(r)e 



-i/3?e fc(r) 



The kernel £> a b(r,r') introduced above is defined as 

% b (ry)=pv- 1 (r,r')8 ab + p 2 g{r)8(r-r') ! 



(10) 



-ee V 2 5(r-r'). 



where v 

Equation ^ carries complete information about the mutual coupling between counte- 
rions and the surface charge disorder. The grand-canonical "free energy" of the partially 
annealed system is then obtained as 



3? = -k B T' In 2f. 



(11) 



The special cases of purely quenched and purely annealed disorder follow from Eq. i ll I l l for 
n — > and n = 1, respectively (see Appendix IbTi. 

Note that here the number of replicas, n = T /T' , has a direct physical meaning of tem- 
perature ratio [28 29 35 ,36]. A close examination of Eq. ^ indicates that the partially 
annealed disorder gives rise to quadratic surface terms of the form g(r)^ a (r)^,(r). It may 
thus lead to renormalization of the mean surface charge (Appendix [C]l as can be seen most 
clearly by looking at the mean-field equations which we shall derive next. 



3 Mean-field limit 



The mean-field or Poisson-Boltzmann (PB) equation 16,7,91 (which becomes exact in the 
limit of small coupling parameters corresponding, for instance, to low counterion valency or 
low surface charge density I11II12I ) follows from the saddle-point equation of the functional 
integral <[Sj as 

ee V% = aqeoQ(T)e- i ^ e ^ r Hi P o{T)+pg{T)^h(r)- (12) 

b 

We shall assume no preferences among different replicas on the saddle-point level, thus 
ty a (r) = (j)(r) for a = 1,...,« (replica symmetry ansatz). In this way we arrive at the PB 
equation for the real-valued mean-field potential 9pb(i*) = i^(r) as 

ee V 2 (ppB(r) = -Xqe Q Q(r) e -^^ -p eff (r), (13) 

where 

p eff (r) = p (r)-j8'g(r)(pp B (r) (14) 

is the effective (renormalized) macroion charge distribution (Appendix |CJ, It is therefore 
seen that in the quenched limit (« or jS' — > 0), the disorder effects completely vanish on the 
mean-field level and the PB theory coincides exactly with that of a non-disordered system 
of bare charge distribution po(r) 1171 . This is however not true for the partially annealed 
disorder (n > 0). 

To proceed with the PB theory, we shall consider the specific case of two parallel charged 
plates located (normal to z axis) at z = — a and z = +a at the separation distance d = 2a (Fig. 
[TJ. We take the mean charge distribution and its variance as 

p (r) = -CTe [S(z + a) + 8(z-a)] (15) 
g(r)=ge 2 [S(z + a) + S(z-a)], (16) 

where g > and without loss of generality we assume that a > (and thus q > 0). Counte- 
rions are assumed to be confined in between the plates (i.e., Q(r) = I for |z| < a and zero 
elsewhere), where Eq. d 1 3b admits the well-known solution [6,7 ,9] 

<PPB(z) = 7^1ncos 2 (/i: Z ) (17) 
pqeo 

with K 2 = 2nl-ftq 2 X to be determined from the electroneutrality condition stipulating that 
the total charge on the two surfaces should be equal to the total charge of the counterions. 
This leads to the equation for K of the form 

Ka tw(Ka) = 1 + ylncos 2 (£a) = — . (18) 

a 

Here p. = \/{2itqino) is the Gouy-Chapman length, £ B — e\j (4neeok B T) the Bjerrum 
length, and 

Orff = o + P' ge (pp B (a) (19) 

the renormalized surface charge density (Eq. dl4b). The latter expression clearly reflects 
the mixed boundary conditions encountered here, resembling the situation in the classical 
charge-regulation problems (5J. The dimensionless parameter 
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gives a measure of the disorder annealing and is obviously proportional to the ratio n = T /T' 
of the counterions and surface disorder temperatures. 

The PB pressure, fpB, acting between the plates is obtained from the standard definition 
jSPpb = w PB (zo) - 5j8eeo(d<Pp B /dz) 2 | :o HEEI for an arbitrary \zq\ < a as 

PPPB -(KM) 2 - (2D 



27id B a 2 



The counterion number density profile between the plates, «pb(z) = A e Pi e o<PPB(z) [6]|7]|9)> 
is obtained as 

nrB^W^M 2 (22) 
27ie B a 2 \cosKzJ ' 

It follows from Eq. d 1 3 b that the mean-field renormalized surface charge density is al- 
ways smaller than the bare value and tends to zero but never changes sign as y increase 
(0 < (7 e ff < c). Therefore, the surfaces are effectively neutralized and the pressure as well 
as the counterion number density profile tend to zero as / increases to infinity (see Fig. [2}. 
This picture relies on the assumption that the number of surface charged units is not fixed 
and can respond to changes of the surface potential. Imposing the constraint that fixes this 
number obviously rules out surface charge renormalization and one observes no effects from 
the disorder annealing in agreement with Ref. 1191 . 

In the limit y — > 0, we recover the non-disordered 1121 or quenched |17 | mean-field 
results with the following asymptotic behavior for the pressure, 

lKi *° 2 \n 2 ii 2 /d 2 rf/M»l- 
For large y » 1, we find that (Ka) 2 ~ (jU/a + y) -1 and thus 

;8^b 



2fx/d yrf/jit<l, 

(24) 



2x£b<J 2 [4^/(^2) y£ // M »i. 
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The small separation expression above is nothing but the ideal-gas osmotic pressure of coun- 
terions that dominates over the energetic contributions. At large separations the pressure is 
found to decay asymptotically as ~ 1/ {yd 2 ). The pressure remains always non-negative and 
the surface-surface interaction is thus always repulsive in the mean-field limit. Choosing 
the non-disordered system as the reference, however, the decrease in the interaction pres- 
sure upon increase of the surface disorder annealing can be interpreted as being due to an 
effective disorder-induced attraction whose asymptotic form could be described by 

1-y 1 

AP PB ~ — ^ (25) 

for large y. This asymptotic form again attests to the fact that the way the disorder acts on 
the mean-field interaction between the two apposed surfaces is via a renormalization of the 
surface charge density. 



4 Strong-coupling (SC) limit 

Next we shall investigate the asymptotic strong-coupling limit which is complementary to 
the mean-field limit and where counterion-induced correlations become dominant. We em- 
ploy the standard strong-coupling scheme reviewed extensively in Refs. II 111121 in order to 
study the partial annealing effects in the SC limit. The so-called asymptotic SC theory is 
obtained from the leading order terms of a non-trivial virial expansion (in powers of the 
fugacity) of the partition function i.e. 

2? = 2f Q +k2f l + ff(X 2 ). (26) 

It becomes exact in the limit of large coupling parameters corresponding, for instance, to 
high counterion valency or high surface charge density II I III2I . 

The zeroth-order (no counterion) term, 3?q, and the first-order (single counterion) term, 
T\ , follow from Eq. l[8} as 

^=/(n5r)^°> (27) 

Si = £ / dRI2(R) / ( fl ^) e --*>-*/»«>*W (28) 

where 

= W /drdr> fl (r)^(r,r')fe(r')+^£ /drp (r)^(r). (29) 

1 a.b J a ■' 

We thus need to calculate both these terms for an arbitrary number of replicas, n. In doing 
so, we shall make use of some mathematical relations that we briefly discuss below. 



4. 1 Mathematical preliminaries 

First, it turns out that the most convenient way to carry out the calculations is to replace 
the long-range Coulomb interaction v(x) = 1/(4jt££o|x|) with the exponentially screened 
Yukawa interaction 

v s (x) = (30) 
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by introducing a finite screening length fc . In the end, we shall take the limit K — > 0. Note 
that not only is this procedure of technical convenience but it is also of physical relevance for 
the present problem. It corresponds to adding a background salt to the system, leading to the 
screened Coulomb interaction between charged units. It is to be noted however that the salt 
effects are taken into account in this way only on the linear Debye-Hiickel level. (Further 
study of the role of added salt in the SC limit is presented elsewhere [27].) Generalization 
of Eqs. ©-(10) in the presence of Yukawa interaction is immediate as v 1 (r,r') is simply 
replaced by 

v; 1 (r,r') = ee (-V 2 + ,c 2 )5(r-r'). (31) 
Second, in calculating 3?q and 3f\ one needs to evaluate the determinant and the inverse 
of the block-matrix £> a i,(r,r'). These calculations are straightforward and may be carried 
out most easily by employing properties of block-matrices and the operator algebra defined 
over the Hilbert space {|r)}. We shall use the compact notation by defining the operators 
(r|v,s|r') = Vj(r,r'), (r|g|r') = g(r) 5(r,r') (for the screened Coulomb interaction and dis- 
order variance), and (r|D a /,|r') = S> a i,(r,r') via Eq. d 10b . where the latter is defined as an 
element of the n x n operator matrix 

D = J Se®v7'+j8 2 u®g (32) 

with e a b = 8 a b and = 1. Also, we shall use (r|po) = po(r) and the well-known notation 
Jr.r'Po(r)v s (r,r')po(r') = (p |v s |p ), etc. 

One can prove the following identities for D 

detD = (det/3v7 1 )"det(l+n J Sgv s ), (33) 

P£(r|(6- 1 mi / ) = (r\v~ s (i+npg v - s )- 1 \r'), (34) 

a 

P (r | (D 1 ) aa | r) = (r | v s [l - pgv s (1 + rajSf v s ) _1 ] | r) , (35) 

which will be used in what follows. Note that the last two quantities do not depend on the 
replica indices a and b. 



4.2 Virial terms % and 

Going back to the virial term one can perform the Gaussian integral in Eq. d27| l to obtain 

% = Coe -ilndet l 6C J .-ilndetD-^ l 6 2 L,./,<Pol(D- 1 U|Po). ( 36 ) 

Using Eqs. ((33) and 434| >. Sq is completely determined. This term represents the interaction 
free energy of macroion charges in the absence of counterions. 
Next, the Gaussian integral in Eq. d28b can be evaluated as 

Jfj =n% JdRn{R)e- p < R \ (37) 

where u(R) is nothing but the single-counterion interaction energy with the macroion charges 
and reads ^ 

u(R)=l3qeo^( Po \(t>- l ) ab \R)+l5 { -^(R\(D- l ) hb \R). (38) 

a z - 

This is fully determined by virtue of Eqs. d34b and d35 b - We have thus derived the general 
form of both virial terms in Eq. d26l l as a function of n. We now proceed to the explicit 
evaluation of the two virial terms for the system of two apposed charged planar surfaces as 
defined by the mean charge distribution and disorder variance H51 and d 16b (see Fig.[T}. 
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4.3 Small-n expansion 



To proceed further we focus on the small-n limit of the above expressions. Our chief goal 
here is to examine the stability of the system upon small annealing perturbations of a 
quenched charge distribution. The annealing effects on this level are therefore expected to 
be additive in the free energy of the system. 

By expanding Eq. d36b for small n we arrive at 



, n/3 
■Qexpj - — 



Tr(gv s ) + (po|v s |p ) 



+ 



(np) 



-Tr({gvJ 2 ) + (po|v s gv 5 |p ) 



(39) 



to the lowest orders in n. But since we are interested in the inter-plate interaction, we shall 
need to determine only the separation-dependent terms. 

It easily follows that the Tr(gv s ) term in the above equation does not depend on the 
inter-surface distance d = 2a and the Tr({gv s } 2 ) term may be calculated straightforwardly 
(and up to an irrelevant additive term) as 

J S 2 Tr({gv ! } 2 ) = -5(4^ 2 g 2 )Ei(-2^), (40) 

where S is the total area of each surface and Ei(x) = J^dte'/t is the exponential-integral 
function. The remaining two terms in the expression for 3?o are obtained as 

J S(po|v. s |po) = 25(a 2 4)^)(l+e^ £/ ), (41) 
/3 2 (p |v s gv s |po> =2S(ga 2 £l) (^\ x+e -*Jf. (42) 

The small-n expansion for 2?\ leads to evaluation of g~^"( R ) ~ e -0"o( R ) ^ -\- n [iu\ (R)] , 
where we have expanded w(R) = «o(R) — nu\ (R) + 0(n 2 ) to the lowest order in n, and thus 
the following two terms 

«o(R) = qe {po\%\R) -J3^|^(R|v^ s |R>, (43) 
«i(R) =p> (po|m|R} - J S 2 ^^(R|v. s gv s |v J |R). (44) 



We have discarded a self-energy term (q , eo) 2 (RK' ! s|R} in Eq. d43b whose only effect is to 
rescale the fugacity, A 1171 . We can then use the explicit expressions 

/ko(Po|v s |R> = ~{o£ B ) (jfj {e-«\°- R t\ + (45) 
^ 2 eo (po|v,|v,|R> = -(ogf B ) (^)\l +e -^)(e-^ +e-^\). (46) 
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The two remaining expressions in Eqs. d43 b and ( 144b are obtained as 



J3 2 4(R|v,|v,|R) = -{Ingll) 



Ei(-2Jc|a-.K z |)+Ei(-2Jc|a+R z |) 



(47) 



p 3 4(R\0 s gv s gv s \R) = {2Kg 2 £l) ( ^ 



e 



-2x\a-R z \ _j_ e -2nr|fl+R-| J r 2 e ~ lKd \ + 



K 

+2K\a-R z \Ei(-2K\a-R z \)+2K\a+R z \Ei(-2K\a+R z \) + 



+4Kd Ei(-2Kd) 



(48) 



We shall need only the small K results, which read (after discarding irrelevant constants) 

J3 2 4(R|v,|v s |R) ~ -{2ngf B )\n(a 2 -R 2 ), (49) 



P 3 4(R\W.§9,\R) ~ (2ng 2 4) ( f ) (e- 2 ^+e- 2 ^ + 2^). 



(50) 



This completes the explicit evaluation of both virial terms. We now proceed to the evaluation 
of the interaction free energy, i.e., the part of the free energy that explicitly depends on the 
separation between the two surfaces. 



4.4 SC free energy 



By using Eqs. d36b- d38b . one can evaluate the free energy of a partially annealed system 
from J? sc = -k B T'\n2? (Eq. CD), where J" = J° + ^^i in the SC limit (Eq. {26)). 
The fugacity can be fixed by the number of counterions N upon transforming to canonical 
ensemble via 

nN = X d -^f, (51) 

whereby we obtain 

N 



JdRi2(R)e-/ i «(R)- (52) 

The canonical SC free energy is obtained via the Legendre transform, = J^ sc + 
Nk B T\nX,as 

^f = _ln^,_ 1 r fl(R)e - W R) (53) 
N Nn J 

supplemented by the electroneutrality condition, which for the two-plate system reads Nq = 
2a e ffS, and again stipulates that the total charge on the surfaces equals the charge of the 
interposed counterions. Note that here CT e ff is the effective surface charge density that has to 
be determined self-consistently within the SC theory. It turns out however that, in the SC 
limit, there is no charge renormalization due to partial annealing on the leading order in n, 
and thus CT e ff = CT up to corrections of the order &{n 2 ) (see Appendix [C}- This behavior is 
in stark contrast with the one found on the mean-field level in Section[3] 

The above expressions d3 1 b— d5 3 b together with Eqs. d36b- d38b are applicable to any gen- 
eral system of fixed macroions with (Gaussian) disordered charge distribution and arbitrary 
degree of annealing n. For small annealing perturbations, we have 

^ = -^- m /dR^(R) e -^^ (54) 
N Nn J 1 ' H JdRI2(R)e-/ i "o(R) v v 
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where is given by Eq. 

In what follows, we shall focus again on the two-plate model system (Eqs. d!5l l and d!6l l) 
and make use of the explicit expressions d40b-d50b in order to calculate the SC free energy 
of this system. We then take the limit of small inverse screening length, K — > 0, as noted in 
Section RTTI We thus find that the SC free energy of this system adopts a simple form as 

- /quenched + ~ (Si? qifago) d , 

/ 2y \ d 

— /quenched ~t~ I I 5 (55) 

where the first term on the right hand side, /quenched, is nothing but the quenched (n = 0) 
rescaled free energy 1171 

/quenched = J~ + (X - *) hd, (56) 

and the second term in Eq. d55l l represents the leading order correction from partial annealing 
of the disorder. Note that this term is linear in n as expected and scales with the inverse 
screening length jc . 

The quenched free energy fl56l l is expressed in terms of the dimensionless disorder cou- 
pling parameter % = 2nq 2 l\g, which gives a measure of the disorder-induced coupling via 
the surface charge variance g. This is to be compared with the so-called electrostatic cou- 
pling parameter E = litq^ll^o defined originally for non-disordered systems |12], which 
measures the strength of counterion-induced correlations and, in the present context, de- 
pends on the mean charge density a. The information about the disorder annealing enters 
only via y = ng/(qo) as defined previously in Eq. d20l >. These dimensionless parameters 
may be used to determine the phase behavior of a partially annealed system, with the disor- 
der effects being in general quantified by % and y. 



4.5 Instability and collapse transition 



The quenched free energy d56l l comprises the standard non-disordered SC contributions, i.e., 
the counterion-mediated attraction, d/2fi, and the repulsion, — lnd, due to the confinement 
entropy of counterions between the two surfaces 11111121 . But it also includes a long-range 
additive logarithmic term, that is ;£lnd, stemming from the disorder variance which is at- 
tractive and renormalizes the repulsive entropic term. This peculiar form of the quenched 
disorder contribution leads to the previously predicted 1171 continuous collapse transition at 
the threshold % c = 1 between a stable bound state of the two surfaces and a collapsed state 
where the surfaces are in contact. In other words, the optimal surface-surface separation, d*, 
behaves as 

d (2(1-*) %<l, 

H (57) 

If on the other hand the disorder is partially annealed (w > 0), we see from Eq. d55l l that 
the annealing generates a linear attractive term in the free energy, that is ~ yd/ K. It adds up 
with and enhances the counterion-mediated attraction term (first term in Eq. fl56l l) exhibiting 
thus a complementary effect when compared to the quenched disorder contribution. The 
threshold of the collapse transition % c = 1 remains intact up to small corrections of the order 
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0{n) 1371 and the partially annealed bound-state separation is obtained in the limit of small 
inverse screening length, K", as 



which is always smaller than the quenched value, reflecting again the enhanced surface- 
surface attraction in the case of partially annealed disorder. 

Note that the optimal separation in the partially annealed case, Eq. d58l l. tends to zero, 
d„ — > 0, as the inverse screening length tends to zero K — > 0. In other words, the system 
goes into a collapsed state regardless of other parameters exhibiting thus a global attractive 
instability. Hence, a quenched system of macroion charges and counterions is not stable 
with respect to annealing perturbations of the macroion charge distribution in the absence 
of screening effects. The stability may be achieved by adding a finite amount of added salt. 
Unfortunately, the effects of added salt have not yet been properly analyzed in the context 
of the SC theory (see Ref. [27 1 for a recent attempt) and it is presently difficult to go beyond 
the linear description adopted here for the salt screening effects. 



5 Conclusion and discussion 

In this work we have analyzed the effects of partially annealed disorder in the distribution 
of macromolecular (macroion) charges on the interaction between two such macromolec- 
ular surfaces across a solution containing mobile neutralizing counterions. Recent experi- 
ments on decorated mica surfaces [2] covered with a random mosaic of positive and nega- 
tive charged domains (stemming from the adsorption of cationic surfactants) clearly point 
to the existence of a strong attractive surface-surface interaction which, upon formation of 
such domains, collapses the system into a compact state with the surfaces being in contact. 
This behavior resembles the transition to a primary minimum within the DLVO theory (6] 
|7) although the attraction mechanism here is strictly non-DLVO |2] producing attractive 
electrostatic forces that are up to a few orders of magnitude larger than the universal van- 
der-Waals forces (§1 as incorporated in the DLVO theory. Moreover, the emergence of an 
attractive instability or even a collapse transition is not predicted within the standard the- 
ories of electrostatic interactions between charged macromolecular surfaces (even between 
surfaces bearing opposite charges) pointing thus to the possible role of the surface charge 
disorder in the aforementioned collapse transition. 

The AFM investigation of the surface texture of mica surfaces (2) can be performed 
only before the measurements of the inter-surface forces and can not be monitored while 
the two surfaces are brought closer together. One can thus not be sure whether the surface 
distribution of the charged domains along the surfaces changes on approach of the surfaces 
or not. This leads in general to three possible scenarios: 

- the surface charge disorder is completely set by the method of preparation of the surfaces 
and does not change on approach of the surfaces (quenched disorder), 

- the surface charge disorder responds to the changes in the separation of the surfaces just 
as fast as the mobile charged species (such as counterions) in the solution between the 
surfaces (annealed disorder), 




(58) 
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- and lastly, the surface charge disorder does respond to the changes of the separation 
but with a much larger time-scale than the mobile charged species between the surfaces 
(partially annealed disorder). 

As the experiment alluded to above gives only the inter-surface interaction as a function 
of the separation, one is thus lead to investigate the fingerprint of these different scenarios 
on the behavior of the interaction. Since the effects of annealed | 5, 19,20,21 ,22,23] and 
quenched 1 17, 18, 19] disorder have already been investigated theoretically, we follow up 
on these analyses by investigating the changes in the inter-surface forces wrought by the 
intermediate case, where it is assumed that there is a clear separation of relaxation time 
scales between the dynamics of the "external" surface charges (T s ) and the mobile charges 
floating in the solvent between the surfaces (T c j), i.e. 

T s > T ci . (59) 

Thus, the fluctuations in the intervening Coulomb fluid have time enough to relax to their 
local equilibrium state for each configuration of the surface charge disorder, which shows a 
much slower relaxation with changes in the surface separation. The origin of the different 
relaxation times for the surface and bulk dynamics could be manyfold: the finite mobility 
and mixing of charged units on the macromolecular surfaces, such as appears in lipid bilay- 
ers with embedded charged proteins |4|, conformational rearrangement of strongly charged 
DNA chains such as appears in DNA microarrays [1 ] and charge regulation of contact sur- 
faces bearing weak acidic groups in aqueous solutions |5], but should be in any case very 
system specific. This separation of time scales is closely related to the so-called "adiabatic 
elimination" of fast degrees of freedom [31 ,32,33], which has been used frequently in the 
literature for systems with widely different time scales [ 28 , 29 , 34 , 35 , 36 1 (see also Appendix 

In this work we thus investigate the "interaction fingerprint" of the partially annealed 
disorder on two apposed planar surfaces upon their approach. The analysis presented above 
points to the fact that partial annealing of the surface charges invariably leads to additional 
attractive interactions between the surfaces and may even result in a global attractive insta- 
bility in the system. The nature of these attractions is quite different, however, if counterions 
are strongly or weakly coupled to charged surfaces in the sense of Netz II 111121 . The mag- 
nitude of this coupling essentially depends on the valency of counterions, q, the magnitude 
of the surface charge density, Geo, and the Bjerrum length £g = e\/{AiZEEQk B T) (incorpo- 
rating the medium temperature and dielectric constant), and is measured by the electrostatic 
coupling parameter 

E = 2nq 3 i\a. (60) 

For weakly coupled counterions, i.e. specifically in the mean-field limit (E — > 0), the 
partially annealed disorder leads to smaller mean-field repulsions due to a renormalized 
(reduced) value of the surface charge density. The attraction in this case can thus be inferred 
only from a diminished repulsion with respect to the case of a non-disordered surface charge 
distribution. For strongly coupled counterions, i.e., on the SC level (S — > °°), we derive 
explicitly an additional inter-surface attraction stemming from the partial annealing of the 
surface charge distribution. Note that the qualitative difference between the mean-field and 
the SC results goes back to the strong electrostatic correlation effects that are included in 
the SC theory but excluded on the mean-field level. The reason for this additional attraction 
(or reduction of the inter-surface repulsion in the mean-field limit) compared to the purely 
quenched case is that any rearrangement of the macroion charges, such as is assumed in 
partial annealing, inevitably leads to configurations of lower (free) energy which shows up 
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as an effective attraction. This means that every disordered charged system will be unstable 
against annealing. In both the mean-field as well as the strong-coupling limit, the effect of 
partial annealing is quantified by a single partial annealing parameter 

qO 

which depends on the temperature ratio n = T /T'. This is in addition to the disorder coupling 
parameter 

X = 27tq 2 £ 2 B g, (62) 

which measures the spread of the charge disorder distribution via the disorder variance g. 

The parameter y and the effective surface temperature T' may be estimated experimen- 
tally from the mobility (TJ) and the diffusion (D s ) coefficients of the surface charges and by 
applying Einstein's relation k^T' = D s /r s (Appendix I A. 11 1. The temperature ratio quantify- 
ing the effects of partial annealing can thus be cast into an equivalent form 

«=(£ B r)if. (63) 

How does one differentiate the effects of quenched and partially annealed disorder? 
Comparing the results derived in the quenched case 11711181 with those obtained in this 
work, one realizes that in the former case the effects of the disorder are limited to the strong- 
coupling limit, while in the partially annealed case they persist also in the mean-field limit. 
This is because in the mean-field limit the surface charge density is renormalized in the 
presence of partially annealed disorder. This effect is always absent in the quenched disorder 
case (Appendix let. 

Moreover, we find that in the SC regime and in the limit of low screening K — > 0, a 
charged system that may be stable in the presence of quenched charge disorder could become 
globally unstable and collapse upon partial annealing of the disorder. One should also note 
that for small non-zero screening and in the stable phase (x < 1), the system adopts a much 
smaller surface-surface distance in the partially annealed case than in the quenched case as 
the ratio between the optimal distances in these two cases is given by ( 1 + Ay/ Kfi)~ l . 

Another indicative feature of partially annealed disorder is the curious property that 
the surfaces interact electrostatically even if the mean surface charge density, a, goes to 
zero, that is when the surfaces become net electroneutral! It may be seen most clearly from 
the Tr({gv,} 2 ) term in Eq. < f40b , which in the limit k — » leads to a logarithmic attractive 
contribution in the free energy as ~ ting 1 (\S lnd stemming only from the variance, g, of the 
surface charge disorder (note however that when a > this is a subleading term compared 
to the additive 1/fC term considered in Eq. d55l l) 1 37 1. This effect disappears in the quenched 
limit (77 — > 0) unless either the dielectric discontinuities at bounding surfaces or the presence 
of salt in between the bounding surfaces are taken into account [18]. It is thus safe to say 
that the effects of partially annealed disorder are in general stronger and ubiquitous and may 
be qualitative in both mean-field and strong-coupling limits. 

Nevertheless, though the above-mentioned interaction fingerprints would help in assess- 
ing the importance of disorder or even its presence in the charge distribution on macromolec- 
ular surfaces, the interactions by themselves would not be enough to make this conclusion 
with a reasonable degree of certainty. Unfortunately, more detailed experiments where one 
would concurrently measure interactions as well as surface charge distributions would still 
be essential. 
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Our analysis of the effects of macromolecular charge disorder on their interactions sup- 
plements the recently acquired new wisdom of bio-colloidal interactions 1 1 1 , 12, 13 14|, as 
opposed to its classical formulation (6]|7), in quite an illuminating way. Whereas on the 
DLVO or mean-field level 16,7,9] one can formulate the salient features of macromolec- 
ular Coulomb interactions with the folk wisdom that opposites attract and likes repel, the 
strong-coupling paradigm IHII12l fT3 ,14] suggests that likes attract too if the system is highly 
charged. To this we would add, in view of our previous work I17II18I and the work described 
above, that if the surface charge distribution is disordered, the system may become unstable 
and collapse due to attractive disorder-induced forces and that even neutral macromolecu- 
lar surfaces can interact via electrostatic interactions. This is especially true if the disorder 
distribution is partially annealed as we set to prove in this work. 
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A Dynamics of disorder 

A.l Adiabatic elimination of fast variables 

The effective partition function for a partially annealed system may be obtained via an adiabatic elimination 
of fast variables discussed widely in the literature [31 32 33 |. Here we shall briefly discuss the extension of 
this method for a charged system with a slow disordered component. As described in the text we shall focus 
on a non-equilibrium system of mobile fast-relaxing counterions (located at positions for i = 1, . . . ,N, 
in a solvent of fixed temperature 7") and a disordered ensemble of macroion surface charges (located at 
positions x, with i = 1, . . . ,N S and bearing charges q s eo) which undergo a much slower dynamics. We shall 
be interested only in the stationary-state properties of the system. We adopt the Langevin equation for the 
dynamics of counterions and surface charges as 

dH 

§i = -r^ + Ba, (64) 

x ; = -r s — + S , (65) 
ox, 

respectively, where the Hamiltonian H is written in terms of Coulomb pair interactions 

H=Wv(x„x J )+1v(x,^ J ) + Wv^,^^ (66) 

Z U i.) U 



and c j and fl s are white noise terms with zero average and the two-point correlation function fld(f) • c i ('') = 
6D cl S(t—t') and 8 s (f) ■ B s (t') =6D S S (( — ('). We impose Einstein's relations between the diffusion constants 
and the mobilities of the particle species, i.e., D cl = k^Tr cl and D s = k^T'r s that guarantee a stationary state 
may be achieved at long times. (We have neglected hydrodynamic interactions in Eqs. (64} and (64) for 
simplicity as they do not affect the stationary-state results.) 

The central assumption for the time scales of the counterions (1/id) and the surface charges (l/r s ) 
translates into 

r s /^.i«i. (67) 

Our goal is to eliminate the the fast (counterionic) variables that change at short time scales and obtain an 
effective dynamical equation for the slow (surface charge) variables assuming that the effective temperature 
of the latter, T', may in general be different from that of the former, T [28, 29 34 35 36]. This may be 
done readily given the condition (67) above and by developing a perturbative expansion (much in the spirit 
of the Born-Oppenheimer approximation) [31 -32, 33! in terms of the small parameter r s /r c i. The zeroth- 
order results are quite intuitive: one ends up with a Langevin equation for x, similar to Eq. (65) where the 
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microscopic Hamiltonian H is replaced with an effective one, which is obtained by "pre-averaging" over the 
counterionic degrees of freedom; namely, 

x , = -r s ^£ + e S! (68) 

where // c ff = — k&T~inZa and Z c ; = /n;d§, exp(— fiH). This effective Hamiltonian is of course nothing 
but the equilibrium free energy of counterions in the presence of a fixed realization of surface charges; its 
derivative acts as a driving force pushing the surface charge dynamics to reach a stationary state at long 
times with the Boltzmann-type probability weight exp(— p"// c ff) as dictated by Eq. (68). Hence, the system- 
although far from equilibrium as a whole-may still be described in terms of an "effective" partition function, 
Z, once a stationary state is achieved, i.e. 

(Zci)", (69) 
where n = /3' /p\ 



dxj 



kin 



dx, 



A. 2 Dynamical density functional theory 

It is possible to derive the «-replica partition function (8) directly from the expression (69) and by using a Leg- 
endre transformation to grand-canonical ensemble described by the partition functions S = s Z and 

3f c i = Y.n ^*Z c j, where A s and A denote the fugacities of the surface charges and counterions, respectively. 
We shall however take a detour over the so-called dynamical density functional theory that first transforms 
the effective Langevin equation (68) into a dynamical equation for the surface charge density field p(r,f). 
This may be done straightforwardly following the approach proposed for classical fluids in Ref. [38]. One 
thus obtains 

^p(r,t) = -V-J p (r,f) + !j(r,0, (70) 
i conserved density field 1 39 40], where the current density reads 



assumine a c 



Jp(T,t) = -D s V P (r,t)-(q s eor s ) P (r,t)V{Y(T,t))T- 



(71) 



The first term on the right hand side represents the current due to the diffusion of surface charges and the 
second term represents the contribution from the mean effective electrostatic force acting on them; here 
(yf(r,f))r = — SlnSci[p]/ fiSp(r) is the electrostatic potential at the mucroion surface averaged over the 
counterionic degrees of freedom. The last term is a Gaussian white noise with T](r,t) = and r) (r,i) r)(r',t') = 
2V- V'[D s q s e p(r,t)8(r-T , )]8(t-t'). Combining Eqs. (To) and {7j], one arrives at 



dt' 



p(r,f) = V • D s Vp + (q s e r s )p V(y}r + 1? (r,( 



( 9s e r s )p(r,f)V 



Sp(v,t) 



SW[p] 



+ l(r,f), 



(72) 
(73) 



where the effective "potential", ^[p], follows as 



^/ dl " 



p In 



1 



lnr a [p] 



(74) 



with po being a constant. Using a standard mapping of Eq. (73) to an equivalent Fokker-Planck equation, one 
can easily show that the stationary-state fluctuations of p(r,f) follow a Boltzmann-type weight exp(— fi'W) 
and hence, the effective "partition function" 



- /dr[pln(p/po)-p] 



SUM')] 



(75) 



In order to proceed with an analytical investigation, we assume that the surface density fluctuations 
around po are small and expand the entropic term (first term in Eq. 174V ) up to the second order, which gives 
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p ln(p/po) — p ~ (p — po) 2 /2po- One can generalize this result by assuming that different realizations of the 
surface charge disorder follow a Gaussian probability distribution as 

.^[p]=Cexp(^-i|dr g - 1 (r)[p(r)-po(r)] 2 j, (76) 

where C is a normalization factor and g(r) is the disorder variance. The quadratic ansatz is completely general 
since it corresponds to an expansion in small deviations from pair) where g (r) can be interpreted as a disorder 
"compressibility". One thus recovers Eqs. and (4j in the text. 



A.3 Non-conserved disorder dynamics 

The dynamical model discussed above is based on a conserved model for the dynamics of the disorder, where 
the current density follows from the gradient of an effective chemical potential, i.e.,Jp = — (g s eo-Ts) p V/i c ff [p] 
with 

as seen from Eq. {73). Alternatively, one could construct other dynamical models that include a non-conserved 
density field for the surface charge disorder and still give rise to the same effective partition function in the 
stationary-state limit [ 39 40j . One may thus propose a phenomenological equation of the form 

j t P(T,t) = -fs(/irff[p] - Mo) +Tj(r,0, (78) 

where r s is an effective mobility, is a constant (which can be set to zero) and rj is a Gaussian white noise 
with T){r,t) = and Tj(r,f)Tj(rV') = 2D S <5(r — r') S(t — t'). Assuming D s = kBT'f s one can show that the 
partition function \75\ is reproduced as the stationary-state solution. This model may be more appropriate to 
study dynamics of disorder in macromolecules with ionizable surface groups. 



B Partially annealed free energy 

The "free energy" of a partially annealed disordered system follows from Eq. @ as 

& = -k B T' )n 3? = -k B T' In ((%£)). (79) 
Obviously, the purely annealed disorder limit with true thermal equilibrium, T = T', is recovered when n = 1, 

^annealed = In «jT ci )), (80) 

while the purely quenched disorder limit follows by taking the limit n — > (or p" — > 0). One can use {{iKf )) 
r * and expand both sides of this relation to obtain the standard quenched free energy expression 1281291 

^quenched = -fc B r«ln^oi)>- (81) 

All other values of n represent a partially annealed situation. The quenched limit translates into an infi- 
nite temperature for the disorder [28 29 36 | relative to counterions, reflecting the fact that the statistics of 
quenched disorder is unaffected by counterions. It also reflects the zero mobility of surface charges when the 
dynamical interpretation via Einstein's relation is used (assuming a finite surface diffusion coefficient). 



C Effective surface charge density 

The effective surface charge density o^ff in the two-plate system is defined via 



(82) 
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where [[—]] = ^ f&p »\p\ (5°d[p])" (•■•)" thefu11 ensemble average with respect to the partition 3?, Eq. 
g). It may be written in terms of ({■ • ■ )) averages defined after Eq. g) as [[■ • ■ ]] = {(( 3f c i[p}) "(•■•)»/« (3Ta [p]' 
One can easily show from Eqs. g) and j 



IP(r)]]=PoW 



J that 

-«(r)|^=Po(r)-i/3g(r)£[[^(r)]]. 



(83) 



Obviously, the effective or renormalized surface charge is shifted from the bare value, Po(r), by an amount 
proportional to the mean surface potential [[0„(r)]]. Thus the charge renormalization in the present context 
is rendered via a linear surface screening mechanism (as may be seen also from the quadratic disorder terms 
g(r)0 a (r)0j,(r) in the effective Hamiltonian (9))- This can be traced back to the Gaussian approximation for 
the disorder distribution, Eq. (5). 

It follows immediately from Eq. {83) that no charge renormalization occurs in the quenched limit (n — » 0) 
and it may be possible only for a partially annealed disorder (n > 0). Assuming the replica symmetry ansatz, 
one can recover the PB effective charge, Eqs. (141 and )19k in the mean-field limit. In the SC limit, the 
effective charge and the fugacity, A, are determined simultaneously from Eqs. 482t and )83t above combined 
with Eq. j5U or (52) in the text and the electroneutrality condition Nq = 2o" c ff5, i.e. 



A = 



9 /dRfl(R)«-M*) ' 



Ccff 



2Se J 



dr 



Po(r)+g(r 



S In 3T 



Spo(r) 



(84) 
(85) 



By using the SC results for the single-particle interaction, u, and the partition function, 3f, from Sections l4.2l 
and !4.3l and expanding g c ff and A in powers of n as 



o- cff = o" (0 ' +na v 
A = A(°)+«A("+- 

we find that O"' ' = a, l (0) = 2oS/q JdRI2(R)«-' i "o( R ', and 
P 



t(1) 



a 

AW 
A(°) 



2o"5'eo . 
CT 



| drg(r) ( PQ \v s \r)+qe kW f dRi2(R) £ -' ! "»( R )(r|v I |R) 



, JdRX2(R) Hl (R)e-' ! "»( R » 
/dRfl(R)e-P"o( R ) 



(86) 
(87) 



(88) 
(89) 



The A ( 1 ' term and other higher-order terms lead to higher-order contributions in the effective charge and thus 
will not be relevant here. By straightforward calculations and using the results in the text, one can show that 



gW (2n\ 



JdRI2(R)e-' s "o( R )(e- ,c l a -^l +<r I 



JdRX2(R)e-' ! "o( R ) 



(90) 



In the limit K — > 0, the first two terms in the bracket on the right hand side are cancelled by similar contribu- 
tions from the integrals in the last term and one can show that 



g« 
lim = 0. 

k^o a 



(91) 



Therefore, in the SC limit, the surface charge density is not renormalized on the leading order in n when 
the low screening limit is taken. Thus, we can impose the electroneutrality condition using the bare surface 
charge density as 2gS = qN. 
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